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n 
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Hi  Undersize  population  density,  number/f t-cu.f t.  of  mother 
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Abstract  of  Thesis  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the  Degree  of 
Master  of  Science  in  Engineering 

THE  STABILITY  OF  A  MIXED  SUSPENSION  CRYSTALLlZER 
WITH  CLASSIFIED  PRODUCT  WITHDRAWAL 

By 

Ashish  Jayant  Mehta 

March,  1969 

Chairman:   Thomas  M«  Reed 

Major  Department:   Chemical  Engineering 

The  effect  of  size  classification  on  the  dynamic  behavior  of  a 
mixed  suspension  crystallizer  with  classified  product  withdrawal  is 
investigated  by  analyzing  a  simple  model  of  such  a  crystallizing 
unit. 

It  is  observed  that  the  stability  of  the  crystallizer  depends  on 
the  relative  kinetic  rates  of  nucleation  and  growth.   The  results  show 
that  classification  lowers  the  sensitivity  of  the  crystallizer  to  dis- 
turbances and  reduces  the  tendency  toward  cyclic  fluctuations . 
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CHAPTER  I 
INTRODUCTION  TO  CRYSTALLIZATION 

Theoretical  Approach 

In  continuous  crystallization  processes,  a  considerable  importance 
is  attached  to  the  technically  complex  problem  of  controlling  the  crystal 
size  distribution.   The  form  of  the  size  distribution  depends  on  the 
physico-chemical  characteristics  of  the  crystallizing  material  and  on  fhe 
constraints  imposed  upon  the  crystallizing  unit  itself.   Industrially  and 
otherwise,  the  study  and  control  of  crystallization  processes  and,  hence, 
of  crystal  size  distribution  have  been  carried  out  for  several  years  [1]. 
However,  even  up  to  the  present  time,  most  of  the  process  techniques  are 
still  something  of  an  art  and  depend  to  a  considerable  extent  on  the  ex- 
perience and  intuitive  skill  of  the  operator. 

Theoretical  approach  to  industrial  crystallization  has  been  late  in 
coming,  and  papers  written  on  the  matter  deal  mainly  with  the  steady  state 
behavior  of  operating  systems.   In  1956,  Saeman  [2]  derived  equations  for 
the  steady  state  size  distribution  in  a  mixed  suspension.   He  compared  the 
theoretical  distribution  with  that  obtained  from  an  Oslo-Krystal  type  unit 
producing  ammonium  nitrate  crystals  under  conditions  of  hindered  settling, 
and  found  a  fairly  good  correlation  between  the  theoretical  prediction  and 
experimental  results.   A  similar  work  was  done  by  Robinson  and  Roberts  [3] 
who  obtained  the  theoretical  crystal  size  distribution  resulting  from  a 
series  of  well-mixed  units  with  nucleation  in  the  first  stage  only.   A 


comparison  of  the  size  distribution  obtained  from  a  single  unit  with 
that  obtained  from  an  industrial  ammonium  sulfate  producing  crystallizer 
showed  good  quantitative  agreement. 

The  special  feature  that  distinguishes  continuous  crystallization 
processes  from  other  continuous  reactors  is  the  simultaneous  occurrence 
of  nucleation  and  growth.   This  point  was  brought  out  clearly  in  the 
work  of  Saeman  and  in  that  of  Robinson  and  Roberts;  and,  indeed,  it  has 
been  observed  that  the  control  of  nucleation  rate  is  the  key  to  con- 
trolling the  size  distribution  in  a  crystallizer.   Saeman  [2]  has  dis- 
cussed the  problem  of  controlling  the  nucleation  rate  by  removal  of 
crystal  nuclei  in  order  to  control  the  size  distribution  for  the  case 
of  a  crystallizer  with  classified  product  withdrawal. 
Dynamic  Behavior 

Papers  dealing  with  the  dynamic  behavior  of  crystallization 
processes  are  yet  fewer  in  number.   This  is  not  surprising  in  view  of 
the  fact  that  the  processes  of  nucleation  and  growth  in  a  crystal- 
lizing system  are  of  a  complex  nature  and  do  not  lend  themselves  easily 
to  mathematical  description.   In  1962,  Randolph  and  Larson  [4]  derived 
time  dependent  equations  relating  crystal  population  density  to  size  in 
a  mixed  suspension  of  arbitrary  volume.   These  equations  are  based  on 
the  concept  of  conservation  of  particle  population,  and  it  was  pointed 
out  that  a  particulate  system  must  conserve  particles  in  addition  to 
satisfying  laws  describing  the  conservation  of  mass  and  energy.   Ana- 
lytically, the  concept  of  population  balance  is  of  importance  inasmuch 
as  it  provides  a  way  of  characterizing  the  size  distribution  by  taking 
into  account  the  effect  of  nucleation  on  the  size  distribution.   The 
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steady  state  form  of  the  population  balance  equation  was  solved  for  the 
case  of  a  series  of  crystallizers  with  nucleation  in  each  unit.   The 
size  distribution  thus  obtained  was  reduced  to  that  obtained  by  Robinson 
and  Roberts  [3]  when  nucleation  was  assumed  to  occur  in  the  first  stage 
only.   Further,  for  the  case  of  a  single  unit,  the  size  distribution  was 
found  to  be  identical  to  that  derived  by  Saeman  earlier  [2]. 

In  another  study,  Randolph  and  Larson  [5]  solved  the  transient 
population  balance  for  the  case  of  an  ammonium  sulfate  producing  mixed 
suspension  crystallizer .   The  transient  response  of  the  size  distri- 
bution to  upsets  in  nuclei  dissolving  rate  and  in  production  rate  was 
studied  by  simulating  the  system  equations  on  an  analog  computer.   The 
results  were  interesting  inasmuch  as  it  was  observed  that  the  size 
distribution  undergoes  long  term  cyclic  fluctuations  due  to  pertur- 
bations affecting  the  system.   In  an  earlier  study.  Miller  and  Saeman  [6] 
reported  similar  observations  in  an  industrial  ammonium  nitrate  producing 
crystallizer.   In  the  field  of  continuous  polymerization  and  in  fermen- 
tation processes,  such  observations  have  also  been  made  [7,  8]. 
Problem  of  Stability 

In  view  of  the  observed  dynamic  characteristics  of  crystallizing 
systems  involving  long  term  cyclic  fluctuations  in  particle  size  distri- 
bution, the  question  of  the  stability  of  such  systems  becomes  one  of 
interest  inasmuch  as  it  has  a  bearing  on  the  problem  of  controlling  the 
size  distribution.   Randolph  and  Larson  [5]  have  investigated  the  sta- 
bility limits  of  a  mixed  suspension  crystallizer  with  mixed  product 
withdrawal.   These  authors  point  out  the  importance  of  the  relative 
kinetic  rates  of  nucleation  and  growth  in  defining  the  stability  limits 
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4 
of- the- system.  •  The  stability  criterion  was  expressed  In  terms  of  the 

ratio  of  nucleation  rate  to  growth  rate  which,  as  the  authors  point 
out,  would  be  transgressed  in  the  event  of  a  discontinuous  jump  in 
the  nucleation  rate.   Such  a  situation  would  arise  for  instance  when 
the  supersaturation  becomes  so  high  that  sudden  mass  nucleation  occurs 
as  the  me tas table  solubility  boundary  is  crossed. 
The  Classified  Product  Withdrawal  Crystallizer 

■   Many  crystallizers  are  operated  with  some  kind  of  a  classification 
device,  either  in  the  crystallizer  bed  itself  or  in  the  outlet  of  the 
crystallizer.   The  types  of  classifiers  in  general  use  are  either  those 
that  segregate  particles  by  gravity  separation  or  those  that  use  centri- 
fugal action.   An  Oslo-Krystal  type  unit  may  be  cited  as  an  example  of 
the  first  type  wherein  only  the  larger  particles  that  settle  down  at  the 
bottom  of  the  crystallizer  bed  are  withdrawn.   Devices  using  the  elutri- 
ative  principle  belong  to  the  same  class.   On  the  other  hand,  devices 
such  as- hydroclones,  which  make  use  of  the  centrifugal  principle,  fall 
into-  the  second  category.   Theoretically,  an  ideal  classifying  device 
would  cause  a  perfect  separation  of  oversize  and  undersize  particles. 
However,  in  most  of  the  known  classifying  devices  in  which  a  certain 
amount  of- mother  liquor  is  withdrawn  to  maintain  overall  material  balance, 
a  proportionate  amount  of  undersize  particles  cannot  be  prevented  from 
being  withdrawn  with  the  mother  liquor.   Further,  a  factor  that  prevents 
the  product  from  being  perfectly  uniform  is  that  product  withdrawal  takes 
a  finite  amount  of  time,  and,  therefore,  the  oversize  particles  exhibit 
a  certain  degree  of  size  distribution.   An  ideal  classifier  would  require 
an  infinite  recycle  rate  for  an  instantaneous  withdrawal  of  oversize 
crystal  product. 
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The  steady  state  aspects  of  a  crystailizer  with  classified  product 
withdrawal  have  been  studied  previously  [2,  9,  10].  Miller  and  Saeman  [6] 
observed  that  classification  tends  to  produce  crystals  of  a  poor  quality. 
They  found  that  it  was  impossible  to  operate  a  completely  classified 
ammonium- nitrate  producing  unit  continuously,  but  that  it  was  necessary 
to  withdraw  the  product  periodically  in-  a  cyclic  manner. 

Sherwin  et  al.  [11]  have  analyzed  the  dynamic  behavior  and  the 
■  stability  limits  of  a  crystailizer  with  ideal  classifying  action.   These 
authors  found  that,  for  the  classified  case >  the  stability  criterion  is 
considerably  relaxed  in  relation  to  the  case  of  mixed  product  withdrawal, 
or  in  other  words,  as  was  pointed  out,  size  classification  introduces   • 
greater  instability  Into  the  system.  However,  the  validity  of  their  re- 
sults is  questionable  inasmuch  as,  for  reasons  given  earlier,  the  concept 
of  an  ideally  classified  crystailizer  is  not  realistic  when  applied  to 
practice.  The  mathematical  description  of  such  a  system,  wherein  there 
is  no  withdrawal  of  undersize  particles  and  the  withdrawal  rate  of  over- 
size particles  is  infinite,  is  not  suitable  for  a  realistic  evaluation  of 
the  actual  situation.  The  present  work  attempts  to  analyze  the  stability 
problem  of  a  classified  product  withdrawal  crystailizer  with  a  more  real- 
istic model  wherein  both  the  undersize  and  the  oversize  particles  are 
removed  at  a  finite  rate,  and  there  is.ai,  preferential  withdrawal  of 
oversize  particles.  A  linearized  stability  analysis  of  the  system  e- 
quations  is  made,  and  the  analytical  approach  follows  closely  that  for 
a  crystailizer  with  mixed  product  withdrawal  as  given  by  Randolph  and 
Larson  [5]. 

The  theoretical  background,  including  the  concept  of  population 
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conservation,  is  developed  in  Chapter  II.   The  stability  problem  of  the 
mixed  withdrawal  crystallizer  is  also  discussed,  and»in  Chapter  III,  the 
classified  withdrawal  case  is  analyzed.   In  Chapter  IV,  the  linearized 
system  equations  are  investigated  for  the  stability  criterion,  and  the 
results  are  discussed  subsequently  in  Chapter  V. 


CHAPTER  II 

THE  POPULATION  BALANCE  AND  THE  MIXED  SUSPENSION 
MIXED  PRODUCT  WITHDRAWAL  CRYSTALLIZER 


Conservation  of  Particulate  Entities 

In  Chapter  I,  it  was  mentioned  that  the  primary  factor  controlling 
the  particle  size  distribution  in  continuous  crystallization  processes 
is  the  kinetics  of  nucleation  and  growth.   Now,  inasmuch  as  the  product 
crystal  size  distribution  is  also  dependent  on  the  residence  time  proba- 
bilities of  particles  within  the  system,  a  theoretical  development  at- 
tempting to  describe  the  size  distribution  in  various  systems  must  take 
into  account  the  interaction  between  the  residence  time  distribution  of 
particles  and  the  kinetics  of  nucleation  and  growth. 

Now,  a  particulate  system  must  conserve  numbers  in  addition  to  con- 
serving mass,  momentum,  and  energy  if  the  size  distribution  is  to  be 
fully  defined;  and,  therefore,  equations  of  the  conservation  of  the  last 
three  quantities  are  not  sufficient  for  formulating  the  problem  of 
particle  size  distribution.   In  general,  one  can  speak  of  the  distri- 
bution of  particulate  entities  over  any  associated  property.   This 
property  can  be,  for  instance,  the  size  of  a  particle,  or  age,  or  any 
quantity  that  varies  as  a  function  of  time.   A  population  balance  of 
particulate  entities  derived  with  sufficient  generality  can  be  used  co 
describe  the  property  distribution  in  various  constrained  systems 
[12,  13,  14].   Such  a  general  derivation  of  a  population  balance  due 
to  Randolph  [15]  is  given  in  the  sequel. 
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Consider  the  distribution  of  particulate  entities  that  are  continu- 
ously distributed  through  a  region  U  of  n-dimensional  space 


\(Xi,  X2,  ...,  x^) 
and  through  a  region  V  of  m-dimensional  property  space 


A  point  number  density  in  the  region  R  =  U  +  V  can,  therefore,  be  de- 
fined as 

n  =  n(xi,  X2,  ...,  x^,  li^i,  11^2,  ••.,  ^^) 

Next,  the  velocities  of  propagation  of  the  entities  along  the  1 

spatial  axis  and  along  the  1   property  axis  may  be  expressed  as 

8x,         9i(;, 
1       _   i 

\.    ~   9t   '  '^^.    ~   9t 
1  1 

Further,  consider  the  generation  and  dissipation  functions  operating 
throughout  the  region  R 
Generation: 

G  =  G(xi,  X2 x^,  \l)i,    ^2>    •••>  ^^>    t) 

Dissipation: 

D  =  D(xi,   X2,    ...,   x^,    i>i,   1|J2,    ...,    ^^,    t) 

which,  of  course,  are  functions  of  time.   Now,  a  particulate  balance 
may  be  written  as 

{Accumulation}  =  {Input}  -  {Output}  (2-1) 

If  a  Lagrangian  approach  is  taken,  that  is,  if  it  is  assumed  that  the 


region  R  translates  with  local  velocities  v   ,  the  particulate  balance 

^i 
Equation  (2-1)  becomes 

1^   TndR  =   TcG  -  D)  dR  (2-2) 

R        R 
where 


dR  =  dxj  dx2  • .  .  dx  dipi  d4/2  •••  d^ 


or  equivalently, 


d_ 
dt 


/ndR  =  ^//      III  I      ^'^\   '"    '^'^^   ^^n   •'  '^^1 

R  Xi  X2  .  .  .  x^  ij;i  ijj2  •  •  •  'l^jj^  (2-3) 


Equation  (2-3)  may  be  differentiated  by  repeated  application  of  the 
Leibnitz  rule  (See  Appendix  I)  to  give 


n  m 


r  {|ii+y  ^(v^,n)  +y  ^(v^,n)  ^jj_c}  dR  =  o      (2-4) 


R      '  i=l  i=l 

which  further  reduces  to 

n  m 


In + y  ^(\.^^  +  y  ^^^.""^ + D  -  G  =  0         (2-5) 

i=l  i=l 

inasmuch  as  R  is  an  arbitrary  region.   Equation  (2-5)  is  a  population 
balance  for  a  particulate  system  but  is  of  a  form  that  is  too  general 
for  application  to  the  size  distribution  problem  in  a  crystallizer , 
However,  before  Equation  (2-5)  is  constrained  to  represent  population 
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density  distribution  in  a  crystaliizer,  it  is  essential  to  discuss 
the  role  of  nucleation  and  growth  in  a  crystallizing  system. 
Nucleations  and  Growth 

The  driving  force  in  the  process  of  crystallization  is  the  super- 
saturation  of  the  solute  material  in  suspension.   Further,  inasmuch  as 
the  nucleation  rate  and  the  growth  rate  both  are  functions  of  super- 
saturation,  one  would  expect  a  functional  relationship  of  the  following 
.type,  expressing  the  dependence  of  nucleation  rate  and  growth  rate  on 
supersaturation . 

Growth  rate: 

'  • ■        dN° 

^=f(s)  (2-6) 

Nucleation  rate: 

r  =  g(s)  (2-7) 

where  N  is  the  total  number  of  nuclei  per  unit  suspension  volume.   The 
functions  f (s)  and  g(s)  must,  of  course,  be  determined  experimentally. 
Kinetic  Growth  Models 

Miller  and  Saeman  [6]  found  a  linear  function  of  the  form 

r  =  kis  (2-8) 

describing  the  dependence  of  growth  rate  on  supersaturation.   The  work 
of  Saeman  [2]  and  of  Van  Hook  [16]  has  supported  this  kind  of  relation- 
ship. 
Kinetic  Nucleation  Models 

Volmer  found  an  Arrhenius  type  kinetic  law  [17]  describing  nu- 
cleation behavior.   His  model  is  based  on  thermodynamic  considerations 
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and  may  be  expressed  as 


dN°   , 

dr  =  '^2  ^ 


X,n(s  +  C  )1 
s 


(2-9) 


where  k2  and  k3  are  constants,  and  C  is  the  solubility  concentration. 

s 

Another  nucleation  model,  which  is  a  simplified  representation  of 
Equation  (2-9),  is  Mier's  metastable  model  [17] 

f^=k,(s-sy  (2-10) 

where  k^  is  a  constant,  and  s  is  the  critical  supersaturation  below 

m 

which  the  solution  is  metastable.   There  is  some  experimental  evidence 
supporting  this  model  [4]. 

At  this  point,  it  should  be  mentioned  that  the  process  of  nu- 
cleation in  a  mixed  suspension  crystallizer  is  a  rather  complex  phe- 
nomenon involving  both  homogeneous  and  heterogeneous  nucleation,  the 
latter  process  being  the  formation  of  nuclei  on  a  crystal  surface. 
Volmer's  model  was  derived  for  homogeneous  nucleation  only,  and, 
therefore,  its  applicability  to  a  mixed  suspension  crystallizer  is 
questionable.   In  a  heterogeneous  nucleation  process,  the  nucleation 
rate  should  depend  both  on  concentration  and  total  crystal  surface.   In 
the  presence  of  a  large  number  of  crystals,  the  dependence  on  surface 
area  is  relatively  less  important  than  the  dependence  on  supersaturation. 
Rumford  and  Bain  [18]  have  investigated  the  nucleation  of  sodium  chloride 
in  a  bed  with  large  crystals.   These  authors  found  that  the  dependence 
of  nucleation  rate  on  supersaturation  is  very  similar  to  that  of  homo- 
geneous nucleation.   For  such  a  system,  Volmer's  model  would  well 
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express  the  process  of  nucleation  in  a  crystallizing  system.   The  ad- 
vantage of  using  Equation  (2-9)  is  that  it  is  smooth  over  the  entire 
concentration  range. 

Nucleation  experiments  on  an  ammonium  sulfate  system  have  indi- 
cated that  the  metastable  supersaturation  level  in  such  a  system  is 
almost  nonexistent  [4].   In  this  event,  Mier's  nucleation  model  would 
reduce  to 


iN°  _  ,,  P 


dt 


=  k^sP  (2-11) 


where  the  value  of  p  has  been  given  as  4  for  the  ammonium  sulfate 
system.   Sherwin  et  al .  [19]  have  offered  criticism  for  the  use  of 
Equation  (2-11)  in  dynamic  studies  on  the  grounds  that  in  a  large 
number- of  systems  metastable  supersaturation  cannot  be  neglected.   It 
will,  however,  be  shown  in  the  sequel  that  this  problem  does  not  arise 
in  determining  the  interactive  relationship  between  nucleation  and 
growth  since  in  a  well-mixed  crystallizer  both  the  nucleation  rate  and 
growth  rate  depend  on  the  same  level  of  supersaturation. 
Nucleation  and  Growth-  Rate  Interaction:   The  Power  Law 

•    The  nucleation  rate  can  be  related  to  the  nuclei  density  n   (which 
is  actually  the  population  density  of  particles  of  "zero"  size)  as 
follows 

dN°   3L  ,3N°,     o  n   17^ 

where  r  =  3L/9t  is  the  crystal  growth  rate  expressed  as  the  rate  of 
change  of  the  characteristic  size  L  of  the  particle;  the  characteristic 
size  being  any  one  spatial  dimension  of  a  collection  of  geometrically 
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similar  particles  that  can  be  used  to  characterize  the  size  of  the 

particles . 

Combining- Equations  (2-8)  and-  (2-11)  results  in  the  power  law 

where  k^  is  another  constant.  In  Equation  (2-13) ,  inasmuch  as  the 
nucleation  rate  is  expressed  as  a  function  of  the  growth  rate,  the 
level  of  supersaturation  on  which  the  former  quantity  depends  does 
not  enter  into  the  expression. 

■    The  advantage  of  using  the-  nucleation-growth  power  model  over 
that  of  Volmer  is  clearly  seen  from  the  form  of  Equation  (2-13) 
which  is  a-  mathematical  simplification  of  Equation  (2-9) .   The  ex- 
ponent- p  in  Equation  (2-9)  is  actually  the  sensitivity  ratio  of  nu- 
cleation rate  to  growth  rate  when  both  are  expressed  as  functions 
of ■ supersaturation.   As  was  mentioned  in  Chapter  I,  this  ratio  is  im- 
portant in  defining  the  stability  limits  of  a  mixed  suspension  crystal- 

lizer. 

Substituting  Equation  (2-12)  into  Equation  (2-13)  results  in  the 

expression 

n°  =  kj^r^  (2-14) 

which  is  another  form  of  the  power  law  and  where  i  =  p  -  1.   Equation 
(2-14)  expresses  the  nuclei  density  in  terms  of  the  growth  rate  and  is 
a  convenient  boundary  condition  for  solving  the  population  balance 
equation  as  will  be  shown  later,  and  it  is  known  to  hold  over  a  wide 
range- of  concentrations  for  many  inorganic  systems  [20]. 


?T^^^ 
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The  Mixed  Suspensipn-  Mixed'  Product  Withdrawal  Crystallizer 

Prior  to  Investigating  a  crystallizer  with  classified  product  re- 
moval, the  case  of  a  crystallizer  with  mixed  product  withdrawal  must 
be  reviewed  since  the  problem  of-  the  stability  of  the  former  is  ana- 
lyzed along  lines  similar  to  the  latter. 

The  general  population  balance  Equation  (2-5)  may  now  be  reduced 
further  to  suit  the  description  of  a  mixed  suspension  mixed  product 
withdrawal  crystallizer. 

In  a  crystallizer  of  arbitrary  volume  V,  the  processes  of  particle 
nucleation  and  growth  can  be  characterized  by  three  orthogonal  spatial  • 
coordinates  x,  y,  and  z,  and  a  size  coordinate  L  defining  the  particle 
size.   Hence,  in  this  case,  Equation  (2-5)  reduces  to 

■   1^  +  V.  (vn)  +  l^^'"''^  +  D  -  G  =  0  (2-15) 

where  n  is  the  population  density  in  a  mixed  suspension  crystallizer. 
V  is  the  vector  velocity  of  the  particles  in  the  spatial  region,  and 
r  =  91/ 9 t  is  the  growth  rate. 

If  it  is  further  assumed  that  the  generation  and  dissipation 
terms  are  absent.  Equation  (2-15)  reduces  to 

If +V.(^)  +|^(^^)  =0  (2-16) 

This  equation  can  now  be  Integrated  over  the  volume  V  of  the  crystal- 
lizer, giving 

{In  +  l-Cj^n)  ^  ^  _^   r  ^.  (-^)  ^^  ^  Q  (2-17) 

9t   oL  J 

V 
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Then,  by  the  divergence  theorem,  the  volume  integral  in  Equation  (2-17) 
can  be  converted  to  a  surface  integral 


I     V-(tn)  dV  =   / 


V("n)  dV  =   /  vndS  '  (2-18) 

V  S 

where  the  velocity  v  is  orthogonal  to  the  surface  S.   The  following 
two  contributions  to  this  integral  are  important. 

a.  The  term  due  to  the  change  in  the  volume  of  the  system, 
which  is 

n^  (2-19) 

dt 

b.  The  integral  over  an  arbitrary  nxmiber  of  inflows  and  out- 
flows, which  yields  a  term-  involving  the  flow  rates 


-  J  ".-i 


(2-20) 

1=1 


where  by  convention  inflows  are  considered  positive.  Equation  (2-17) 
thus  becomesf 


n 


9n  +  i_(rn)_y  ^^^+^Aiml=o  (2-21) 

at   3L       /  ^   1  1       dt 


i-1 

The  following  constraints  may  now  be  applied  without  much  loss  of 
generality. 

1.  The  feed  is  clear,  that  is,  there  are  no  nuclei  in  the  feed. 

2.  The  crystallizer  volume  is  constant. 
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Under  these  conditions,  Equation  (2-21)  reduces  to 

3U.  +  i-^^i^)  +  ^  =  0  (2-22) 

at   3L       V 

where  Q  is  the  volumetric  withdrawal  rate. 

In  order  to  uncouple  the  second  term  in  Equation  (2-22) ,  it  is 
necessary  to  assume  a  certain  functional  dependence  of  the  growth 
rate  on  the  size.   Sherwin  et_  al .  [19]  have  assumed  a  growth  rate 
model  of  the  form 

r  =  1  +  aL  (2-23) 

where  a  is  a  constant.   Canning  and  Randolph  [21]  have  used  the  same 
growth  model  to  derive  expressions  for  the  size  distribution  in  a 
mixed  suspension  crystallizer  under  various  constraints. 
Steady  State  Size  Distribution 

Under  conditions  of  a  steady  state  and  assuming  that  McCabe's  AL 
law  [24]  holds,  that  is,  the  growth  rate  is  independent  of  size,  Equation 
(2-22)  can  be  integrated  to  give 

n  =  n°  e"  V  (2-24) 

o    o  '^ 

where  the  subscript  o  refers  to  the  steady  state  condition.   Equation 
(2-24)  suggests  that  in  a  mixed  suspension  crystallizer  with  mixed 
product  withdrawal  the  population  density  decays  exponentially  with 
size.   In  Figure  1,  the  logarithm  of  the  population  density  is  plotted 
as  a  function  of  size.   The  slope  of  the  straight  line  is  -  QL/V,  and 
the  intercept  on  the  y-axis  is,  of  course,  the  nuclei  density  n  . 
Several  systems  are  known  to  exhibit  this  kind  of  behavior  [6,  22]. 
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FIGURE  I :    POPULATION     DENSITY     PLOT    FOR    A    MIXED    SUS- 
PENSION   MIXED    PRODUCT    WITHDRAWAL  CRYSTALLIZER 
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Stability  Analysis 

The  form  of  Equation  (2-22)  suggests  that  even  under  conditions 
of  steady  state  energy  and  material  balance  transients  in  the  size 
distribution  are  expected.   Randolph  and  Larson  [5]  investigated  the 
problem  of  the  stability  of  size  distribution  in  a  mixed  suspension 
mixed  product  withdrawal  crystallizer .   The  linearized  moment  trans- 
forms of  the  system  equations,  obtained  from  the  population  balance 
Equation  (2-22)  in  conjunction  with  the  initial  condition  Equation  (2-24) 
and  the  boundary  condition  Equation  (2-14) ,  were  solved  on  an  analog 
computer.   Step  changes  in  production  rate  and  in  nuclei  dissolving  rate 
were  intrdduced  into  the  system,  and  the  resulting  disturbances  in  total 
suspension  population,  crystal  length,  and  area  were  plotted  as  functions 
of  time.   The  plots  showed  wave-like  cyclic  fluctuations  in  these  quanti- 
ties propagating  with  time,  and  the  periods  of  the  disturbances  were 
found  to  be  three  or  four  times  the  order  of  magnitude  of  residence  times 
of  particles  in  industrial  units.   Transients  have  been  observed  in 
commercial  crystallizers  [6];  and,  in  view  of  this  observation,  the 
question  of  the  stability  of  an  operating  crystallizer  becomes  one  of 
Interest.  The  same  authors  analyzed  the  system  equations  of  a  mixed 
product  withdrawal  crystallizer  by  application  of  the  Routh-Hurwitz 
stability  criterion  and  found  that  the  system  exhibits  instability  when 
p  >  21,  that  is,  when  the  nucleation  rate  varies  with  a  power  of  growth 
rate  greater  than  twenty-one. 

Sherwin  et  al.  [19]  have  also  investigated  the  same  problem.   Their 

treatment  of  the  situation  is  more  general  since  they  have  considered 
the  effects  of  seeded  feed  and  of  nuclei  dissolving  on  the  system 
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stability.   These  authors  have  expressed  the  stability  criterion  in 
terms  of  a  ratio  b/g  which  they  have  defined  as 


d   ,dN°, , 
^dT  ^—^  ^o  ^o 

^/§  "   dr    dN° 

©   (77-) 
ds  o  dt  o 


(2-25) 


Substituting  Equation  (2-12)  into  Equation  (2-25)  gives 

V  pr  P-1  (^)  r^ 
b/g  =  2L.2 ds_g_o  ^  (2-26) 

ds  o  N  o 

For  a  clear  feed,  it  was  found  that  the  system  becomes  unstable  when 
b/g  >  21,  an  observation  which  is  in  agreement  with  that  of  Randolph 
and  Larson.   A  seeded  feed  was  found  to  improve  the  stability  in  certain 
cases.   An  interesting  observation  in  connection  with  the  effect  of 
liquid  voidage  was  that  change  in  the  liquid  void  fraction  does  not  have 
a  significant  bearing  on  the  stability  criterion.   In  a  qualitative 
sense,  this  observation  appears  to  be  correct.   As  it  will  be  shown  in 
Chapter  V,  a  fluctuation  in  the  voidage  causes  a  perturbation  in  the 
rate  of  nucleation  in  the  opposite  sense,  such  that  the  two  effects  tend 
to  nullify  each  other,  so  that  the  dynamic  behavior  of  the  system  is  not 
altered  significantly.   Another  observation  made  was  that  the  system  sta- 
bility increases  significantly  if  the  nuclei  have  a  finite  size  and  that 
increasing  the  nuclei  size  enlarges  the  region  of  stable  operation. 

The  important  point  brought  out  in  these  studies  is  that  the  rela- 
tive kinetic  order  of  nucleation  rate  to  growth  rate  plays  a  significant 
role  in  determining  the  stability  limits  in  a  mixed  product  withdrawal 
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crystallizer;  and  that, for  unstable  operation,  the  criterion  of  p  >  21 
establishes  the  stability  limit  of  the  system  wherein  the  feed  is  clear, 
and  the  nuclei  have  a  size  of  "zero"  dimension. 


CHAPTER  III 

THE  MIXED  SUSPENSION  CRYSTALLIZER  WITH 
CLASSIFIED  PRODUCT  WITHDRAWAL 


Description  of  the  Model 

The  concept  of  population  conservation,  developed  in  Chapter  II, 
may  be  applied  to  a  mixed  suspension  crystallizer  with  classified 
product  withdrawal.   Any  model  describing  such  a  system  should  conform 
to  simple  mathematical  formulation  and,  at  the  same  time,  must  be 
flexible  enough  to  fit  a  realistic  situation.   Consider  a  mixed  sus- 
pension crystallizer  with  a  classifying  device  as  shown  in  Figure  2, 
wherein  particles  below  the  classification  size  L   (undersize)  are 
withdrawn  at  a  rate  Qi,  and  particles  above  the  classification  size 
(oversize)  are  removed  at  a  rate  Q2  such  that 

-     •      Qi  =  aQ2 

where  a  is  the  ratio  of  the  two  flow  rates,  and  Qi  is  also  the  volumetric 
withdrawal  rate  of  the  clear  liquor  for  maintaining  material  balance 
within  the  system.   The  volumetric  feed  rate  is  Qq,  and  the  crystal- 
lizer has  a  constant  volume  V.   The  classifying  device  C  with  a  recycle 
is  shown  as  an  external  component  of  the  crystallizer;  however,  size 
classification  may  occur  Internally  as  in  an  Oslo-Krystal  type  unit. 
In  either  case,  the  basic  mechanism  of  size  classification  remains  un- 
altered and,  thus,  lends  itself  to  the  same  mathematical  description. 
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FIGURE  2 :    THE     MIXED    SUSPENSION     CRYSTALLIZER     WITH 
CLASSIFIED     PRODUCT      V/ITHDRAWAL 
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It  is  further  asstimed  that  the- crystallizer  is  operating  under  conditions 
of  steady  energy  throughput;  that  is,  the  process  is  isothermal,  and  is 
under  a  steady  state  mass  balance  with  constant  material  densities. 
Population  Balance 

Two  population  balances  are  needed  to  describe  the  particle  popu- 
lation distribution  within  such  a  system,  one  for  the  undersize  and  one 
for  the  oversize.   These  are 


9ni     9ni   Qinj 

+  r  -r-  +  -TT-   =0  L  in  [0,  L  ]  (3-1) 


at     9L     V  ^    '     c- 


and 


8n2     9n2   Q2'^2 

where  ni  and  n2  are  the  population  densities  of  the  undersize  and  over- 
size particles,  respectively,  and  where  the  growth  rate  is  assumed  to 
be  independent  of  size. 

Further,  let  r  =  V/Q2  be  the  average  residence  time  of  the  particles 
based  on  the  drawdown  of  oversize  particles.   Equations  (3-1)  and  (3-2) 
then  become 


9ni     8ni   an^ 

+  r  TT-  +  ^r-  =  0  L  in  [0,  L  ]  (3-3) 


9t      9L     r  ^    '      C 


and 


3n2     9n2   n2 

Steady  State  Population  Density 

At  this  point,  it  is  of  interest  to  analyze  the  steady  state 
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behavior  of  Equations  (3-3)  and  (3-4).   At  steady  state,  these  equations 
become 


dn^, 
dL~ 


an 


iO 


r  r 
o  o 


L  in  [0,  L  ] 
'  c 


(3-5) 


and 


dn2o     n2o 


dL 


r  r 
o  o 


L  in  [L  ,  "] 


(3-6) 


which  may  be  integrated  using  the  following  boundary  conditions 


B.C.  1  at  L  =  0,  nio  =  ^lo 


Thus, 


B.C.  2  at  L  =  L  ,  nxo  =  ^20 


-  aL 


r  r 
o    00 
r^lo  -   ^lo  e 


L  (l-a)-L 
c 

r  r 
o      00 
n2o  =  n2o  e 


(3-7) 


(3-8) 


Several  cases  of  these  expressions  are  of  interest. 
Case  1  (g  <  1) 

In  this  case,  as  shown  in  Figure  3,  the  drawdown  rate  of  oversize 
particles  is  greater  than  that  of  undersize  particles.   In  other  words, 
there  is  size  classification,  and  the  slope  of  Equation  (3-8)  is  greater 
than  the  slope  of  Equation  (3-7).   In  Case  a  =  0,  there  would  be  no 
withdrawal  of  undersize  particles  so  that  Equation  (3-7)  would  be  repre- 
sented by  a  horizontal  line.   This  idealized  situation  is  shown  in 
Figure  4. 
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Case  2  (g  =  1) 

This,  of  course,  is  the  same  as  the  case  of  mixed  product  with- 
drawal inasmuch  as  the  rates  of  removal  of  oversize  and  undersize 
particles  are  equal.  Figure  1  depicts  this  case. 

Case  3  (g  >  1) 

This  case  would  arise  if  the  undersize  particles  were  prefer- 
entially removed  from  the  system.   As  is  shown  in  Figure  5,  in  such  an 
event,  the  slope  of  Equation  (3-8)  would  be  less  than  that  of  Equation 
(3-7). 

Case  4  (L  =  0,  L  ^  «) 
c      c 


In  both  of  these  cases,  there  is  really  no  size  classification, 
and,  therefore.  Equations  (3-7)  and  (3-8)  are  reduced  to  those  for  a 
mixed  product  withdrawal. 
Dimensionless  Equations 

The  following  dimensionless  groups  may  be  defined  in  order  to  ex- 
press Equations  (3-3)  and  (3-4)  in  dimensionless  forms 

t       n        L       1  I    r 

6  =  -  ,  y  =  -^  .  X  =  -^Y-   .  and  <!,  -  ^ 

o  o  o  o 

Hence,  Equations  (3-3)  and  (3-4)  become 


9yi    9yi 

+  ^  - —  +  gyi  =  0  X  in  [0,  B]  (3-9) 


89    ^  3x 


and 


3y2     9y2 

- —  +  ())  - —  +  y2  =  0  X  in  [6,  «>]  (3-10) 


86    "^  3x 


where  6  =  L  /r  T  is  the  dimensionless  classification  size, 
c  o  o 
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Moment-  Transformation 

The  n   moment  of  the  population  distribution  is  defined  as  the 


integral 


I  (Population  density)  L  dL,  n  =  0,  1,  2,  ...        (3-11) 


which  for  a  classified  product  withdrawal  crystallizer  may  be  written 
as  the  sum  of  two  integrals 

y  =  /   niL^dL  +  /   n2L"dL,  n  =  0,  1,  2,  ...        (3-12) 


0  L^ 

The  physical  significance  of  the  moments  is  as  follows 

L 

c 


yo  =  /    n^dL  +  /   n2dL  =  The  total  number         (3-13) 
J  J  of   particles  in 


of  particles  in 
L  the  system 


h 

c 


=  j   niLdL  +  / 


yi  =  I   niLdL  +  I   n2LdL  =  The  total  length  of      (3-14) 

all  the  crystals  in 
L  the  system 

L 

c 


/niL^dL  +  k       /      n2L2dL     =  The   total  crystal 
J  area  in  suspensior 


3ion 


^c  (3-15) 


where  k.  is  the  area  shape  factor. 


pk  y3  =  pk  I   n^L^dL  +  pk   /  n2L^dL  =  The  total  mass 
^       ^J  J  of  crystals  in 

suspension 
0  ^c  (3-16) 


29 
where  k  is  the  crystal  volume  shape  factor,  and  p  is  the  density. 

V 

Now  the  weight  distributions  of  the  product  crystals  may  be  written 
as 

Wi  =  pk  niL^  (3-17)  , 

W2  =  pk  n2L3  (3-18) 

so  that  the  mean  crystal  size,  with  respect  to  the  weight  distribution, 
is 

L. 


c 

/  WiLdL  +  /   W2LdL 


0  L 

L  = (3-19) 

L        <» 
c 


/widL  +  j      W2dL 
0        L. 


which  may  be  written  in  terms  of  the  moments  as 


L  (3-20) 

IJ3 


Similarly,  the  variance  of  the  weight  distribution  is 

L 
c 


/   (L-L)2  WidL  +  j    (L-L)2  W2dL 


0 0_ 

a  =  - 


L 
c 


j      WidL  +  r  W2dL 
0  L_ 


(3-21) 
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L          00              L        09  L  " 

c                       c  c 

f    WiL2dL  +  f  W2L2dL  -  2L   {    f    WidL  +  f   WzdL}  +  L 2  {  T  WjdL  +  f  WjdL} 

0           L              0         L^  0  ^ 


L 
c 


r  WidL  +  I       W2dL 
0        L^ 


or  in  terms  of  the  moments 


a (-]  (3-22) 

P  3    'P  3' 

Finally,   the  coefficient  of  variation  may  be  written  as 

(„^,h          U3         '^3' 
C.V.    =  -^ (3-23) 

The  moments  in  Equations  (3-20),  (3-22),  and  (3-23)  can  be  deter- 
mined by  a  sieve  analysis  of  the  product,  and,  conversely,  a  knowledge 
of  these  moments  is  sufficient  to  characterize  the  performance  of  a 
crystallizer.  Further,  as  it  will  be  shown,  the  partial  differential 
Equations  (3-9)  and  (3-10)  can  be  transformed  into  a  set  of  ordinary 
differential  equations  involving  the  moments.   These  moments,  which 
represent  the  total  properties  of  crystal  number,  length,  area,  and  mass 
through  Equations  (3-13),  (3-14),  (3-15),  and  (3-16),  respectively,  de- 
fine the  stability  of  the  system.   Normalizing  Equation  (3-12)  about  the 
steady  state  gives  B  " 


n  normalized 


f  yjx'^dx  +  f  y2x"dx  =  f ^^  +  g^      (3-24) 
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Next,  multiplying  Equation  (3-9)  by  x  and  integrating  over  the  size 
range  x  in  [0,  6]  gives 

6 

/9yi      9yi 
{x"  ^Q-  +  x''  ^^  +  ax^i }  dx  =  0  (3-25) 

0 

which  reduces  to  a  set  of  ordinary  differential  equations  involving 

the  underslze  moments  f  .  This  set  of  equations  is 

n 


dfo 

—  +  4,(Y-y  )  +  afo  =  0 


(n  =  0)      (3-26) 


df 


— r-  +  <I>(3"y  -  nf  ,)  +  af  =0    (n  5^  0) 
de    •'"•'"  I     n-l      n 


(3-27) 


where 


yi 


n°    o 

x=0   n 
o 


and 


yi 


x=e 


Similarly,  Equation  (3-10)  may  be  moment  transformed  to  yield  a 
set  of  ordinary  differential  equations  involving  the  oversize  moments  g^. 


/ 


n  3y2     n  3y2    n 
{x  —  +  (t)x  -^  +  x  y2}  dx  =  0 


(3-28) 


which  gives 


dgo 

^  -  *Y  +  go  -  0 


(n  =  0)      (3-29) 


«ad 
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j^  -  <^(&^y   +  ng^_^)  +  g^  =  0      (n  ?t  0)     (3-30) 

Steady  State  Moments 

Under  conditions  of  a  steady  state,  Equation  (3-24)  becomes 

3  00 

(u  )     ij   J    /yix'^dx  +  /yox'^dx  =  f   +  g       (3-31) 
no  normalized  -J  J  no   *no 

0         3 

The  steady  state  normalized  population  distributions  yi  and  y2  are  the 
dimensionless  forms  of  Equations  (3-7)  and  (3-8).   Substituting  these 
into  Equation  (3-31)  gives 

e        «> 

/   \  r  "Oix  n,   ,   r  B(l-a)-x  n,       ,_  ,-. 

(u)     ,jj    /e    xdx+/e        xdx      (3-32) 
no  normalized  =  J  J 

0  B 

Expanding  the  integrals  on  the  right  of  Equation  (3-32)  gives 

-   -     ,  nl     -a6  f^  ^  ne""^  .  n(n-l)6""^  +  . . . 

f   =  i — 77  -  e  1''=—  H 5 —  +  5 

no     n+1        a     oi'^       a-' 
a 

+  ^  +  ^]},     n  =  0,  1,  2.  ...      (3-33) 

i-.i.^-  a  a 

and 

g   =  e"'*^  [b"  +  ne""-'-  +  n(n-l)  Q^~^   +  ...  +  nig  +  nJ], 
no 

n  =  0,  1,  2,  ...      (3-34) 
so  that  the  first  few  moments  may  be  written  as 

f  oo  =  r  u  "  «  ) 
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fio  =-^U-  e"*^  (l+aS).! 

-)]     (3-35) 


f3o=!r[l-e-«e  (l  +  a6+4^  +  4^'- 


and 


gOo  -  e 

glo  =  e""'^  (6  +  1) 

g2o  =  e'«^  (b2  +26+2) 

g3o  =  e""^  (6^  +  382  +63+6)  (3-36) 

At  this  point,  it  is  interesting  to  note  the  relation  between  the 
solids  residence  time  r  and  the  liquor  residence  time  in  a  crystal- 
lizer  with  classified  product  withdrawal. 

The  average  solids  residence  time  is  defined  as 

crystal  holdup  volume 


s   crystal  volumetric  production  rate 


c 


k^{  r     nioL^dL+r     nzoL^dL} 


0  L 

r     3c  ^    Ml     I     P. 

s  L 

.0 


■—  j;aQ2    j    nioL^dL  +  Q2     T    n2oL3dL}  (3-37) 


L 
c 


which  upon  normalization  of  the  involved  quantities  becomes 
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•^30  +  g3o  I  Y     ff  3o  +  g30 


(^30  ^  630  I   V      I'-iO     •     BiU    I 
1  —  =  ( r  (3-38) 
af 30  +  g3oJ  Q2    »af 3o  +  gsol 

For  the  case  of  mixed  product  withdrawal,  a  is  unity,  hence, 

r  =  r  (3-39) 

s 

or  in  other  words,  the  solids  residence  time  and  the  liquor 

residence  time  are  equal,  which  is  expected. 

Alternatively,  the  production  rate  can  be  written  in  terms  of 

the  growth  rate  and  the  total  surface  area  as  (See  Appendix  II  for 

the  derivation) 

■;;••:.    Production  rate  =  (-r)  r  k^M2oP  (3-40) 

so  that  the  solids  residence  time  may  be  expressed  as 

L 
\{    I      nioL^dL  +  r  n2oL3dL}  p 
0  L^ 

r ^ (3-41) 

s      L  «> 

c 


k^{  f     nioL2dL  +  f   nzoL^dL}  (2^)  p 


0  L^ 

Further,  if  the  shape  factors  k  and  k  are  defined  in  terms 
of  a  unit  cube  then 

k,  =  6k  (3-42) 

A     V 

so  that  Equation  (3-41)  is  reduced  to 

1   1^30  +  830. 

r  =  ■:r-    U 1 (3-43) 

s   3r^  lf2o  +  g2o' 

Now,  for  the  case  of  a  crystallizer  with  an  ideal  classifying  action, 
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Equation  (3-43)  is  further  simplified,  giving 


f30   L, 


r  =  -^  T—  =  7^  (3-44) 

s   3r^  f2o   4r^ 


or 


L  =  4r  r  (3-45) 

c     so 

Thus,  the  classification  size  is  four  times  the  growth  rate  per  solids 
residence  time,  which  was  derived  by  Saeman  [2]. 
Constraint  on  the  Growth  Rate 

The  growth  rate  can  be  expressed  in  terms  of  the  moments  inasmuch 
as  it  is  constrained  by  material  balance.  An  overall  material  conser- 
vation for  the  crystallizer  may  be  expressed  as 

{mass  in}  -  {mass  out}  =  {accumulation}        (3-46) 

Here,  the  accimiulation  term  involves  two  factors,  namely,  the  accumu- 
lation of  mass  in  supersaturation  and  the  accumulation  of  mass  on  the 
crystals.  The  former  factor  is  small  compared  to  the  latter  and, 
therefore,  may  be  neglected.  Thus,  the  material  balance  becomes 

QoCo  -  (ct  +  1)  QiCi  =  (j)  k^yzVp  (3-47) 

and  Equation  (3-47)  reduces  to 

QoCo  -  (a  +  1)  QiCi  =  (|)  k^n°  (r^F^)  3  (fj  +  82)  Vp    (3-48) 


or 


Ko 

r  =  T— (3-49) 

f2  +  g2 


where 
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QqCo  -  (cx  +  1)  QiCi 

c.n°  (r  r  ] 
A  o   o  o 


Ko  °  -I -T (3-50) 

(7)  kV  (r  r  )3  Vp 


Thus,  the  growth  rate  is  inversely  proportional  to  the  total  crystal 
.  area.         ■ 

The  normalized  growth  rate  becomes 

J.     f20  +  g20 

Expressing  the  Quantities  y  and  v  as  Functions  of  the  Moments 

Now  that  the  growth  rate  is  a  known  function  of  the  moments. 

Equations  (3-26),  (3-27),  (3-29),  and  (3-30)  can  be  used  to  solve  the 
stability  problem.   These  equations  form  a  closed  set  at  n  =  2  inasmuch  as 
the  growth  rate  involves  the  second  moments  of  distribution.   However,  the 
quantities  y  and  y°,  or  the  values  of  the  normalized  population  density  at 
the  classification  size  and  the  nuclei  density,  respectively,  are  still 
unknowns  in  these  equations;  and  their  behavior  as  functions  of  the  moments 
must  be  kno^^n. 
Hong's  approximation 

The  following  approximation  expressing  the  relationship  between  y 
and  the  moments  is  due  to  Hong  and  Larson  [23]. 

Consider  the  integral  . 

3  6    3 

/  ""  aT"  ^'^  "  ""^1    "  /  yi^^  =  6y  -  fo      (3-52) 
0  0    0 

wherein  the  right-hand  side  is  obtained  by  the  method  of  integration  by 
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parts.  Further,  changing  the  limits  of  the  function  on  the  left  gives 

e  Y 


/  ^^^^=  /  ^^yi 


(3-53) 


0                          y° 

so  that 

Y 

/      xdyi   =   8y  -  fo 
0 

y 

(3-54) 


The  mean  value  theorem  may  now  be  used  since  x  is  continuous  in  [0,  g], 
and  Yi   is  bounded.  Thus, 

Y 
gy  -  fo  =  J     xdyi  =  (y  -  y°)  X  (3-55) 


o 

y 


X  may  now  be  approximated  as 


^'f7 


so  that 

fl 
f" 


(3-56) 


By  -  fo  =  (y  -  y°)  T-  <3-57) 


and ,  therefore , 


fo  -  y  fl 

Bfo  -  fl 


(3-58) 


which  expresses  y  in  terms  of  the  zero  and  the  first  moment  of  distri- 
bution. Alternatively,  y  may  be  expressed  as  a  function  of  the  oversize 
moments  inasmuch  as 

Y  =  yi(B)  =  y2(3)  (3-59) 
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and,  therefore,  following  the  same  approach 

By  -  go  =  -  YX  (3-60) 


f     9y2 
J  ""IT'  - 


and  letting  x  =  gl/go  gives 


-g^ 


(3-61) 


^  '  (Bgo  -  gl) 

One  way  to  establish  the  validity  of  Equations  (3-58)  and  (3-61)  is  to 
obtain  the  steady  state  value  of  y   by  substituting  the  corresponding 
values  of  the  undersize  and  oversize  moments  in  these  equations.   In 
both  cases, 

Y  =  e""^  (3-62) 

o 

which  is  correct. 

y  and  the  power  law 

In  order  to  express  the  normalized  nuclei  density  y  in  terms  of 

the  moments,  the  power  law  of  Equation  (2-14)  is  used  since  it  is 

assumed  to  relate  correctly  the  nuclei  density  to  the  growth  rate. 

Normalizing  Equation  (2-14)  gives 

y°  =  4,^  (3-63) 

and  the  value  of  ij)  may  be  substituted  from  Equation  (3-51)  to  give 

1^20  +  g2o»i 


y  =  ^-i^ (3-64) 

I  f2  +  g2  ' 
The  stability  problem  can  now  be  analyzed  since  all  the  quanti- 
ties are  now  known  in  terms  of  the  moments  of  distribution. 


CHAPTER  IV 
STABILITY  ANALYSIS 

Linearization  of  Moment  Equations 

In  Chapter  III,  the  basic  population  balance  equations,  describing 
the  distribution  of  crystal  entities  within  a  mixed  suspension  crystal- 
lizer  with  classified  product  withdrawal,  were  converted  into  a  set  of 
linear  differential  equations  in  terms  of  the  moments  of  distribution. 
It  was  also  mentioned  that  a  set  of  these  equations  involving  up  to  the 
second  moments  close  satisfactorily,  and  that  this  set  is  essentially 
sufficient  for  a  stability  analysis.  This  set  of  equations  is 


^+*(Y^-y°)+afo  =  0  (4-1) 


dfi 

j^+   <t.(6Yo  -  foo)  +  af  1  =  0  (4-2). 


df2 

_  +  ,,,(32^^  _  2fio)  +  af2  =  0  (4-3) 


dgo 

—  -  <^Y^  +  go  =  0  (4-4) 


—  -  (t)(6Yj,  +  goo)  +  gl  =  0  (^-5) 


dg2 

J^  -   4>(62y^  +  2gio)  +  g2  =  0  (4-6) 
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where 


U^  -  y°fi 


^  =  Tfrrir  <*-'> 


and 


_   |f2o  +  820  »i 

•^     »  f2  +  g2  ' 

The  set  of  Equations  (4-1)  through  (4-8)  are  now  linearized  by 
Introducing  linear  perturbations  in  the  system  variables.  A  linearly 
perturbed  quantity  may  be  expressed  as  the  sum  of  the  steady  state 
value  and  the  perturbation.  Thus, 


f  (t)  =  f   +  f  (t)  (4-9) 

n      no    n 


8n<^^  =  ^no  -^  ^n^'^  ^^'^^^ 


Y(t)  =y^  +  yit)  (4-11) 

y°(t)  =  1  +  9°(t)  (A-12) 

(t)(t)  =  1  +  *(t)  (4-13) 

where  the  quantities  with  ~  are  the  perturbations,  and  the  steady  state 
values  of  the  normalized  (about  the  steady  state)  nuclei  density  and 
growth  rate,  of  course,  are  unity.   The  linear  equations  therefore  be- 


come 


_  =  _  Y  +  (1  _  Y^)  <!>  +  y  -  af  0  (4-14) 


where 


and 
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dfi      .  .   ~     ~ 

—  =  -  3y  +  (2f lo  -  ey^)  ^  +  ^0  -   af 1        (4-15) 


df2 

^^  =  -  ^2^  +  (2f lo  -  e^Y^)  <|)  +  2f  1  -  af2      (4-16) 


dgo 

_  =  Y  +  Yq*  -  go  (^-^^) 


^  =  By  +  (BYq  +  goo)  *  +  go  -  gi         (4-18) 


dg2     _  -     ~     . 

^^  =  B^Y  +  (a^Y^  +  2gio)  ()>  +  2gi  -  g2        (4-19) 


4,  =  Ki(f2  +  g2)  (4-20) 

y°  =  i^  (4-21) 

Y   =  K2fo  +  K3y°  +  K^f i  (4-22) 


Ki  = J (4-23) 

^   fzo  +  820 

2foo  -  6y„ 
K2  =  T^ T^  (4-24) 

-  fio 

K3  =  -TT T-  (4-25) 

'       Sfoo  -  fio 


K,  =  ,/°  "  f  (4-26) 

It  is  interesting  to  note  that  all  the  quantities  involved  in  Equations 

(4-14)  through  (4-26)  are  dimensionless,  and,  therefore,  the  stability 

of  the  system  can  be  defined  entirely  in  terms  of  dimensionless  quantities. 
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Equations  (4-14)  through  (4-19)  are  linearly  4ependent  inasmuch  as 
the  moments  fg  and  f^  are  related  to  each  other  through  Equation  (4-22). 
Thus,  the  five  independent  equations  may  be  expressed  in  the  vector  form 


d_ 

de 


fl 

fl 

f2 

f2 

go 

=  [A] 

go 

81 

gl 

g2 

g2 

(4-27) 


where  the  transformation  matrix  [A]  Is  obtained  by  substituting  the 
values  of  y  and  y  from  Equations  (4-21)  and  (4-22)  into  Equations  (4-15) 
through  (4-19). 

Tlie  matrix  elements  are  shown  in  Appendix  III.  Th<^  rank  of  the  non- 
symiaetric  matrix  is  five.  In  the  matrix,  the  (iOeff icient  F5  is  the 
ratio  fio/foo- 
Stability  Analysis  of  the  Matrix 

The  matrix  [A]  has  a  fifth  order  characteristic  equation  of  the  form 


X^  +   Ai+A^  +  A3X3  +  A2X2  +  AiA  +  Aq  =  0 


(4-28) 


in  which  the  coefficients  of  the  terms  are  functions  of  the  dimensionless 
steady  state  parameters. 

A  stability  test  can  be  made  by  applying  the  Routh-Hurwitz  criterion 
to  Equation  (4-28).  However,  inasmuch  as  the  process  of  evaluating  the 
coefficients  of  Equation  (4-28)  is  somewhat  tedious,  an  alternative  method 
consisting  of  evaluating  the  roots  of  the  characteristic  equation,  or  the 
eigenvalues,  was  used  in  the  present  study.  As  will  be  seen  in  the  sequel. 
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an  investigation  of  the  root-locus  plots  obtained  from  the  eigenvalues 
is  sufficient  in  defining  the  stability  limits  of  the  system.   The 
digital  analysis  was  made  on  an  IBM  360/50  computer. 


...::^' 


CHAPTER  V 
RESULTS  OF  THE  ANALYSIS 

Root-Loci 

The  root-locus  plots  obtained  for  various  chosen  values  of  the  para- 
meters a  and  6  are  shown  in  Figures  6  through  20.  Only  the  complex  roots 
are  plotted  since  they  are  important  in  determining  the  stability  limit. 
In  a  stable  system,  all  the  roots  lie  on  the  left-hand  side  of  the  com- 
plex plane,  and  a  crossover  to  the  right-hand  side  implies  that  the 
system  is  exhibiting  instability.   In  the  present  analysis,  only  the  com- 
plex roots  crossover  to  the  right-hand  plane  at  the  stability  limit. 

The  following  cases  are  of  particular  interest: 
Mixed  product  withdrawal  (a  =  1) 

Figure  6  shows  that  when  the  value  of  a  is  unity,  the  observed  be- 
havior is  that  of  a  mixed  product  withdrawal  crystallizer ;  that  is,  the 
system  exhibits  instability  at  p  >  21.   The  value  of  the  dimensionless 
classification  size  g  is,  of  course,  arbitrary  in  this  case. 
Size  classification  (a  <  1,  S  =  D 

Figures  7  through  13  show  plots  for  a   values  ranging  from  0.5  to 

0.0.   In  all  of  these  cases,  the  system  is  found  to  be  more  stable  than 

when  the  withdrawal  is  mixed.   It  is  seen  from  these  figures  that  all 

the  root-loci  are  nearly  identical.   From  this  observation,  it  may  be 

concluded  that  the  system  stability  increases  rapidly  when  a  takes  on  a 

value  that  is  less  than  unity,  that  is,  when  the  oversize  withdrawal 

rate  is  larger  than  the  undersize  withdrawal  rate.   However,  little 
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change  in  the  stability  criterion  occurs  once  a  is  less  than  unity. 
Another  point  of  interest  is  the  case  of  no  undersize  withdrawal  (a  =  0) . 
Figure  13  shows  that  even  for  this  idealized  case,  the  system  continues 
to  exhibit  a  high  degree  of  stability. 
Preferential  removal  of  undersize  particle  (a  >  1,  3=1) 

This  case  would  arise  if  the  classifying  action  was  reversed,  that 
is,  if  the  underflow  was  made  higher  than  the  overflow.  Figure  14  shows 
that,  in  this  event,  the  stability  criterion  is  relaxed,  and  the  cross- 
over occurs  at  P  =  17  for  a   =  1.2. 

Effect  of  varying  the  dimensionless  classification  size 
The  following  cases  were  investigated: 

For  a  =  0.005;   g  =  2,  3,  and  4 
For  a  =  0.001;   0  =  2,  3,  and  4 
The  root- loci  for  these  cases  are  shown  in  Figures  15  through  20. 
As  expected  from  the  previous  results  for  values  of  a  less  than  unity, 
the  stability  criterion  is  found  to  be  independent  of  the  value  of  a, 
but  changes  with  increasing  value  of  the  dimensionless  classification 
size.  A  comparison  of  Figures  15  and  18  with  Figures  11  and  12  shows 
that  for  a  given  value  of  a  the  system  stability  increases  greatly  when 
3  changes  from  a  value  of  1  to  2.  For  higher  values  of  B,  there  is  little 
change  in  the  stability  criterion. 
Effect  of  Suspension  Void 

It  would  be  interesting  to  investigate  the  effect  of  fluctuations 
in  the  suspension  void  on  the  stability  of  the  system.  The  void  fraction 
E  in  the  crystallizer  may  be  expressed  in  terms  of  the  moments  as 

e  =  1  -  k^(f3  +  gs)  (5-1), 
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where 

k  =  k  n°r'*r'+  (5-2) 

M    V  o  o  o 

Further,  the  nucleation-growth  rate  law  of  Equation  (2-14)  must  be  re- 
defined inasmuch  as  the  nuclei  density  is  based  on  the  free  volume. 
Thus, 

•    ■   -       n°  =  ek^r^  (5-3) 

Introducing  the  dimensionless  void  fraction       ,       ,    .  j'\ 


V  = 


^  ■  :  :     .     r       (5-4) 


Equations  (5-1)  and  (5-3)  may  be  made  dimensionless  to  give 

1  -  k^(f3  +  g3)        .    ^       ■      ■' 


V  = 


1  -  k^(f3o  +  g3o) 
and 


(5-5) 


■     -         y°  =  v<^^  (5-6) 

The  linearized  forms  of  Equations  (5-5)  and  (5-6)  are 

V  =  K6(f3  +  gs)  (5-7) 

where 

Si 

^''   1  -  V^3o  +  g3o)  ^^"^^ 

and 

y°  =  V  +  i^  (5-9) 

So  that,  substituting  Equation  (5-7)  into  Equation  (5-9)  gives 

y°  =  i^  -  K6(f3  +  is)  (5-10) 
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If  Equation  (5-10)  were  used  instead  of  Equation  (4-21)  in 
stability  analysis,  a  matrix  of  rank  seven  would  result  inasmuch 
as  the  third  moments  of  distribution  are  now  required  to  close  the 
set  of  moment  equations. 
Discussion 

Figure  21  is  a  block  diagram  depicting  the  interrelationship  be- 
tween the  various  quantities  involved  in  defining  the  size  distribution. 
The  level  of  supersaturation  s  is  determined  by  the  feed  rate  Qo  and 
the  total  crystal  surface  area  A.  The  rate  of  nucleation  and  growth, 
in  turn,  depend  on  the  supersaturation.   Finally,  as  is  seen  on  the 
right,  the  size  distribution  is  defined  by  nucleation  and  growth  through 
the  population  balances.  However,  since  the  total  surface  area  is  a 
function  of  the  size  distribution,  there  is  an  internal  feed  back  in 
the  system  as  shown.  This  feed  back  is  the  cause  of  the  cyclic  insta- 
bility in  the  crystallizer  [4]. 

The  important  point  brought  out  in  the  present  analysis  of  the 
stability  problem  is  that  size  classification  tends  to  stabilize  a 
crystallizer.  The  residence  time  probabilities  of  the  undersize  and 
of  the  oversize  particles  are  important  in  defining  the  stability  limit^ 
of  the  system.   If  the  residence  time  of  the  oversize  particles  is  re- 
duced relative  to  that  of  the  undersize,  the  system  shows  increase 
stability.   On  the  other  hand,  if  the  smaller  particles  are  drawn  down 
at  a  relatively  lower  rate,  the  stability  criterion  is  relaxed. 
Sherwin  et  al .  [11]  with  their  model  of  an  ideal  classifying  device 
having  no  drawdown  of  undersize  particles  and  with  an  infinitely  high 
rate  of  withdrawal  of  oversize  particles  find  that  classification 
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actually  introduces  a  high  degree  of  instability  into  the  system. 
They  give  a  value  of  2.3  for  the  sensitivity  ratio  of  nucleation 
rate  to  growth  rate  as  compared  to  a  value  greater  than  21  found  . 
in  the  present  analysis  as  the  criterion  for  stability.  How- 
ever, for  reasons  given  in  Chapter  I,  the  validity  of  their  work  is 
in  question  inasmuch  as  their  model  does  not  describe  the  function 
of  an  actual  classifying  device  in  a  realistic  manner. 

It  has  been  shown  [4]  that  in  a  single  mixed  suspension  crystal- 
lizer,  the  dominant  crystal  size,  with  respect  to  the  weight  distri- 
bution, occurs  at  3=3.  From  the  present  analysis,  it  appears  that 
the  classified  product  withdrawal  crystallizer  operates  with  a  high 
degree  of  stability  when  size  classification  occurs  near  the  region  of 
dominant  crystal  size. 

The  stability  criterion  defined  in  terms  of  the  sensitivity  ratio 
p  would  be  violated,  as  it  was  mentioned  earlier,  in  the  event  of  a 
discontinuous  jump  in  the  nucleation  rate.  This  would  happen,  as  was 
pointed  out,  at  high  supersaturations  when  a  fluctuation  in  the  feed 
concentration  causes  the  supersaturation  to  cross  the  metastable  limit, 
resulting  in  mass  nucleation.   Industrial  crystallizers  do  operate  near 
this  region  of  mass  nucleation  in  order  to  achieve  a  maximum  growth 
rate  consistent  with  low  nucleation  rate.  Mass  nucleation  is  therefore 
likely  to  occur  when  there  are  upsets  in  the  feed  or  in  the  temperature 
of  the  system.  Miller  and  Saeman  [6]  have  made  such  observations  in  a 
crystallizer  with  a  classified  bed  producing  ammonium  nitrate  crystals. 

Referring  again  to  the  block  diagram  of  Figure  21,  it  is  apparent 
that  a  positive  perturbation  in  the  suspension  void  in  the  system  will 
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cause  an  increase  in  the  level  of  supersaturation  inasmuch  as  the 
suspension  area  is  reduced.  This,  in  turn,  will  increase  the  nu- 
cleation  rate  so  that  the  surface  area  will  tend  to  increase  again. 
Qualitatively,  it  therefore  appears  that  the  inclusion  of  the  effect 
of  suspension  void  in  the  analysis  will  not  have  a  significant  bearing 
on  the  stability  criterion. 

Summary 

A  realistic  model  of  a  mixed  suspension  crystallizer  with  classi- 
fied product  withdrawal  was  analyzed  for  stability  limits.  The  moment 
transforms  of  the  population  balance  equations  for  the  model  were  solved 
for  their  characteristic  roots  on  a  digital  computer  for  a  range  of 
classification  sizes  and  at  different  residence  times  of  undersize  and 
oversize  particles.  The  stability  criterion  was  expressed  in  terms  of 
the  sensitivity  ratio  of  nucleation  rate  to  growth  rate,  and  it  was 
found  that,  in  general,  size  classification  tends  to  stabilize  the 
classified  product  crystallizer.  The  following  points  are  relevent  to 
the  current  discussion: 

1.  The  form  of  the  nucleation-growth  rate  law  used  for  the  analy- 
sis is  of  considerable  importance  in  predicting  the  dynamic  behavior 

of  the  system  inasmuch  as  the  stability  limit  of  the  system  is  defined 
in  terms  of  the  relative  order  of  nucleation  rate  to  growth  rate.   Un- 
fortunately, as  was  pointed  out  earlier,  experimental  information  on 
this  important  interactive  effect  is  scarce,  and  the  power  law  used  in 
the  present  study  is  just  an  approximation. 

2.  The  residence  time  was  arbitrarily  based  on  the  drawdown  rate 
of  oversize  particles.   It  is  likely  that  a  more  definitive  way  of 
parameterizing  the  problem  exists. 
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3.   Hong's  approximation  for  correlating  the  population  density 
at  the  classification  size  with  the  moments  of  distribution  may  not 
correctly  represent  the  actual  dynamic  behavior. 

In  view  of  the  above  considerations,  it  is  clear  that  the  present 
work  needs  further  elaboration.  Also  lacking,  among  other  things,  are 
essential  experimental  results  for  gaining  insight  into  the  dynamic 
behavior  of  the  classified  product  crystallizer .  More  information  in 
this  area  is  necessary  for  a  greater  clarification  of  the  problem. 
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APPENDIX  I 
LEIBNITZ  RUI,E 

Leibnitz's  rule  for  the  differentiation  of  definite  integrals 
is  given  as 

b'(t)  b(t) 

a(t)  a(t)  (lA-1) 

where  the  function  f(x,t)  is  integrated  over  the  variable  limits 
a(t)  to  b(t).. 
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APPENDIX  II 

CRYSTAL  PRODUCTION  RATE 

The  total  mass  in  the  crystallizer  is 


I-c        -  ^c 


=  f    WidL  +  f     W2dL  =  k^p  j      niL^dL  +  k^p   T  nzL^dL     (3A-1) 


M 

0  L  0  L 


c 


Now  at  specified  drawdown  rates  and  at  a  given  classification  size, 
the  total  mass  is  constant.  Therefore, 

0  L^ 


or 


9ni  r  8n2 


/dni  f  0112 

l3  —  dL  +    f      l3  -^  dL  =  0  (3A-3) 

0  L^ 

Substituting  the  values  of ■  the- partial  differential  quantities  within 
the  integrals  from  Equations  (3-3)  and  (3-4)  gives 

.  -  .        L_  " 

3ni        ani  r  9n2       n2 


/dnj        ani  r  0112        "2 

l3{-  r  ^^^ p-}   dL  +      /      l3{-  r  -^  -  — >  dL  =  0    (3A-4) 


L 
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These  integrals  upon  evaluation- by  the  method* of "  integration 
by  parts  give 


a    I    niL^dL  +  /  n2L3dL 


0          L. 
r — : £ (3A-5) 

L  =° 

c 


3r   r  niL^dL  +  /  ngL^dL 


L 
c 


So  that  the  production  rate  may  be  expressed  as 

k 

Production  rate 

i 

0  L- 


TT-  {ct  r  niL3dL  +  r  n2L3dL}  =  (|)  rk^yzP      (3A-6) 


where  the  shape  factors  k  and  k  are  expressed  in  terms  of  a  unit 
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